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1.  INTRODUCTION 


Like  most  engineering  design  problems,  the  design  of  a  weapon-vehicle 
system  for  given  objectives  is  first  subjected  to  an  analysis  of  a  model. 
The  model  should  be  such  that  it  describes  the  actual  system  as  closely  as 
possible  and  be  still  within  the  analysis  capabilities.  The  model  of  a 
weapon-vehicle  system  considered  in  this  report  is,  in  some  respects, 
similar  to  the  models  presented  in  previous  Themis  Reports  [1,2].  However, 
it  differs  from  the  previous  models  in  that  it  includes  the  simulation  of 
some  environmental  conditions.  The  primary  interest  is  to  study  the  role 
of  these  environmental  conditions  in  the  design  problems. 

In  Section  2  of  this  report,  the  weapon-vehicle  mode]  Is  described. 

A  detailed  analysis  and  equations  of  motion  of  the  system  are  presented. 
Then  the  problem  is  formulated  as  a  parametric  optimal  design  problem  in 
Section  3.  A  complete  analysis  of  environmental  parameters  and  a  target 
analysis  are  presented  there.  No  numerical  results  are  presented  in  this 
report.  This  is  due  to  the  fact  that  the  problem  of  parametric  optimiza¬ 
tion  has  not  been  explored  in  any  detail  in  the  past.  Currently,  this 
problem  is  receiving  considerable  attention  at  the  University  of  Iowa.  So, 
any  developments  from  this  effort  will  be  reported  at  a  later  time. 
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2.  DYNAMIC  ANALYSIS  OF  THE  WEAPON-VEHICLE  MODEL 

2.1  Description  of  the  Model 

Figure  1  shows  the  model  of  weapon-vehicle  system  considered  in  this 
report.  The  gun  and  its  flexible  mount  are  attached  to  a  vehicle  supported 
by  a  suspension  system  that  is  simulated  by  a  model  similar  to  the  one  pre¬ 
sented  in  [3,4].  Here,  the  suspension  model  is  a  two-dimensional  analog 
of  the  three-dimensional  model  presented  in  [3j.  The  mount  is  considered 
as  an  elastic  structure  with  a  finite  number  of  elastic  elements.  Further¬ 
more,  the  gun  is  assumed  to  have  a  barrel  moving  back  and  forth  due  to  a 
recoil  mechanism  simulated  by  a  reaction  force  R(t) ,  a  spring  and  a  damper. 

The  masses  m^  and  m^  represent  wheels,  axles  and  associated  mechanism 
of  the  rear  and  the  front  of  the  vehicle,  respectively.  The  equivalent 
spring  and  damper  coefficients  for  the  rear  suspension  are  denoted  by  k^ 
and  c^ ,  respectively,  and  the  spring  k^g  is  included  to  denote  a  spring¬ 
like  property  of  the  tire.  As  will  be  discussed  later,  inclusion  of  the 
variable  force  f^  acting  on  m^  will  explain  the  base  displacement  due  to 
various  road  conditions.  Similar  variables  with  subscript  2  represent 
properties  of  the  front  suspension  of  the  vehicle. 

The  main  body  of  the  vehicle  is  idealized  as  a  rigid  body  A  on  which 
the  weapon  C  +  D  is  connected  through  a  flexible  mount  B.  The  weapon  is 
composed  of  a  rigid  housing  assembly  C,  a  rigid  gun  barrel  D,  and  a  recoil 
mechanism  whose  presence  is  denoted  by  the  reaction  force  R(t)  acting  at 
the  ints  Q  and  Q^,  and  a  spring  and  damper  system  kr  and  cr,  respec¬ 
tively.  For  the  purpose  of  representing  other  types  of  recoil  mechanisms, 


the  reaction  force  R(t)  is  assumed  to  be  a  known  function  of  the  relative 
displacement  of  the  gun  barrel  and  its  housing.  The  flexible  part  B  repre¬ 
sents  a  set  of  structural  elements  such  as  truss,  beam,  plate,  and  shell. 

In  the  present  formulation,  B  is  assumed  to  be  a  beam  element  with  both 
ends  clamped  at  the  points  Bfi  and  Bg  of  bodies  A  and  C. 

Following  assumptions  are  made  in  the  analysis: 

1)  The  weapon  is  automatic  with  relatively  high  number  of  bursts. 

2)  The  aiming  by  the  operator  is  perfect. 

3)  Ever}  motion  is  confined  in  a  plane. 

4)  The  whole  system  is  located  in  an  inertial  frame  R  translating  with 
a  uniform  velocity  in  the  plane  considered.  The  notation  R  should 
not  be  confused  with  the  reaction  force  R(t)  or  its  magnitude 
R(t). 

5)  The  point  masses  ,  m^  and  the  mass  center  of  A  can  move  only 
perpendicular  to  a  flat  surface  fixed  in  R. 

6)  The  flexible  mount  B  is  composed  of  linearly  elastic  structural 
elements . 

7)  No  dry  friction  and  no  mass  variation  exist. 

8)  The  axis  of  the  gun  barrel  D  passes  through  the  mass  center  of  the 
housing  assembly  C. 

9)  The  flexible  mount  B  can  be  subjected  to  a  finite  element  analysis. 
This  is  discussed  in  a  later  section  for  a  general  approach. 

Other  assumptions  will  be  made  whenever  necessary  during  the  analysis. 
With  these  assumptions,  the  model  is  a  holonomic  mechanical  system  having 
seven  degrees  of  freedom,  since  infinite  degrees  of  freedom  of  the  flexi¬ 
ble  mount  B  can  be  assumed  to  be  represented  by  a  finite  number  of  boundary 
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displacements.  The  generalized  coordinates  chosen  are;  1)  the  displace¬ 
ments  6^  and  62  measured  from  the  unstretched  lengths  d^  and  d2  of  the 
springs  and  k2Q»  respectively,  2)  the  rotation  in  R  and  linear 

displacement  6^  of  the  mass  center  of  body  A  measured  from  a  configuration 
represented  by  the  unstretched  springs  k^,  k2^,  k^  and  k2»  3)  the  end 
deflection  6^  and  the  slope  02  of  the  tip  of  the  beam  element  with  the 
displacement  along  the  central  axis  neglected,  and  4)  a  relative  displace¬ 
ment  <5r  of  the  gun  barrel  D  with  respect  to  the  housing  assembly  C.  The 
other  necessary  geometrical  dimensions  are  denoted  by  d^,  0,  and  e^ 

1*1,  2,  ...  The  Inertia  properties  of  the  rigid  bodies  are  assumed  to 
be  known  and  will  be  denoted  by  M  and  I  with  appropriate  subscripts. 

2.2  Analysis  of  Flexible  Mount 

An  exact  dynamic  analysis  of  the  flexible  mount  should  be  based  on  a 
continuous  model,  but  for  the  purpose  of  this  report,  a  finite  element 
model  is  used.  For  static  problems,  the  displacements  in  a  continuous 
st.  'icture  can  be  related  to  a  finite  number  of  displacements  selected  at 
some  arbitrary  points  on  the  structure,  usually  the  boundary  poir>':s  inter¬ 
secting  with  the  surrounding  systems.  For  dynamic  problems,  the  displace¬ 
ments  and  velocities  of  any  point  are  implicitly  related  through  the 
differential  equations  of  motion  to  the  history  of  the  displacements  and 
velocities  of  boundary  points  interacting  with  surroundings.  No  explicit 
expression  can  be  found  for  general  dynamic  loadings.  For  a  simple  anal- 
sis,  therefore,  an  approximate  treatment  is  necessary  to  reduce  the  system 
to  a  finite  dimensional  one. 

Recently,  some  approximate  techniques  have  been  used  in  the  analysis 
of  a  flexible  satellite.  In  Ref.  [5], the  authors  used  synthetic 
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modes  in  the  process  of  truncating  the  modal  coordinates  obtained  from  a 
linearized  system,  assuming  small  motions.  Robe  and  Kane  [6]  considered 
two  symmetric  rigid  bodies  connected  by  a  flexible  beam.  They  neglected 
masses  of  the  connecting  structures,  and  linearized  the  equations  with  a 
small  angle  restriction. 

Another  approach  is  to  assume  that  the  accelerations  are  small  and 
use  static  displacement  distributions  for  the  dynamic  case.  This  approach 
is  conventional  [7],  and  is  suitable  for  matrix  dynamic  structural  analysis. 
In  this  method,  the  element  stiffness  and  mass  matrices  for  each  element 
are  generated  in  a  local  coordinate  system.  These  are  then  transformed  to 
a  datum  coordinate  system,  and  combined  to  form  structure  stiffness  and 
mass  matrices.  To  remove  the  internal  degrees  of  freedom,  the  virtual  work 
principle  is  applied.  For  details,  the  reader  is  referred  to  [7,8].  Let 
the  stiffness  matrix  K  and  the  displacement  vector  u  be  partitioned  as 

K11  K12 

K  -  (2.1) 

_K21  K22 

(2.2) 

where  the  vector  u^  refers  to  the  boundary  displacements  to  be  retained  as 
the  degrees  of  freedom  for  the  dynamic  analysis,  and  Uj  denotes  all  the 
remaining  displacements,  i.e.,  the  internal  degrees  of  freedom.  Then  the 
resulting  condensed  stiffness  and  mass  matrices  corresponding  to  u^  are 
given  by 
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K11  ”  K12  K22  K21 


(2.3) 


c  T 
M  -  L  ML 


(2.4) 


where 


\ 

| 


L 


-«22>'1  "22 


(2.5) 


and  1  is  an  Identity  matrix.  Note  that  these  matrices  are  time  dependent, 
since  the  transformation  matrix  from  the  local  system  to  datum  system 
contains  time  variable.  The  Internal  displacements  u 2  are  related  to  the 
boundary  displacements  u^  by 


u 


2 


K21U1 


(2.6) 


In  the  present  problem,  some  further  approximations  will  be  made  concerning 
the  flexible  beam  element.  In  obtaining  a  strain  energy  expression,  the 
static  configuration  will  be  used  and  for  kinetic  energy,  and  potential 
energy  due  to  gravity,  a  lumped  mass  system  will  be  used. 

2.3  Equations  of  Motion  of  the  System 

The  equations  of  motion  of  the  system  will  be  derived  using  Lagrangian 
formulation.  The  notation  used  is  the  same  as  in  Ref.  [9].  For  conven¬ 
ience  of  easy  modifications  of  the  resulting  equations  in  case  of  some 
adjustments  in  the  model  at  a  later  stage  of  numerical  solutions,  the 
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expressions  of  kinematic  analysis  are  written  in  detail  followed  by 
dynamic  analysis  leading  to  the  equations  of  motion. 

2.3.1  Kinematic  Analysis 

For  the  kinematic  analysis  of  the  mechanical  system,  let  (1,  J,  K) 
be  the  right-handed  orthonormal  vectors  fixed  in  the  inertial  reference 
frame  R,  and  let  the  right-handed  orthonormal  vectors  (i^,  j^,  k^)  and 
(i2»  ^2’  ^2^  ke  fixed  in  the  bodies  A  and  C,  respectively,  as  shown  in 
Fig.  1;  here  k^  ■  k2  •  K.  In  order  to  abbreviate  klnematical  equations, 
let  cosine  and  sine  functions  be  denoted  by  c  and  s,  respectively.  Then, 
for  the  base  vectors, 


1  "  c  °1  H  "  8  al 
J  ■  8  i^  +  c 

i2  "  c(®2  +  V*1  +  s(°2  +  6i^i 

j2  -  "  s(a2  +  91)I1  +  c<a2  +  61)J1  . 

The  angular  velocities  of  the  rigid  bodies  are, 

R-A  ,  - 

u  ■  K  , 

or 

CD  -  (^  +  42)K  , 


(2.7) 


(2.8) 


(2.9) 


(2.10) 


i 

< 

» 

i 
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R-D  R-C 

u  ■  0) 


(2.11) 


Now,  the  position  and  velocity  vectors  of  various  points  of  interest  are 
obtained.  A  superscript  *  on  an  entity  refers  to  the  position  of  mass 


center  of  a  body.  For  Body  A, 

rk*/0  -  (d2  +  63) J  , 

(2.12) 

RvA*  »  63(s  a1  ^  +  c  jj)  . 

(2.13) 

For  a  simplified  analysis  assume  that  the  relative  displacements  of  Body 
B  with  respect  to  Body  A  are  small.  From  the  given  end  deflection  6^  and 
its  end  rotation  (*2,  one  can  write  down  the  position  vector  of  B*  as 
follows  flO] : 


_  B*/fi 

r  8  "  e5  ^1  “  ^hl  ^4  +  h2  °2^1  *  (2.14) 


where 


(2.15) 


(2.16) 


vpmm*** 
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The  angle  of  rotation  at  any  point  from  the  original  axis  Is  obtained  by 
differentiation  of  (hj.  6 4  +  h2  <*2)  with  respect  to  e5-  The  velocities  of 
the  location  e5  is  given  by 

*5®*  .  {$3  8  ax  •  hj  ■  h2  &2  -  a^(e^  +  ej)}^ 

+  {6^  c  +  a^(e3  "  ^l  ^4  ”  ^2  a2^^1  *  (2.17 

For  the  point  Be, 

fi  /0  .  —  .  “  /A  |  Q 

r  e  -  (d2  +  63)J  +  (e3  -  6^^  +  (e4  +  e^jj  ,  (2-18 


®v  6  ■  (63  s  *  i4  -  +  e6^^1 


+  (l3  c  ax  +  a^e.^  “  64^Ji  • 


(2.19) 


Similarly,  for  the  Bodies  C  and  D,  one  cub  obtain  the  following  position 
and  velocity  vectors: 

-C*/°  „  -Be/0  _  {e ^  c(a2  +  6l)  +  eg  s(c2  +  6^}^ 

-  {e?  s(a2  +  ox)  -  eg  c(a2  +  ej))?!  »  (2*20) 
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+  Mp  8(65  +  e15  -  e9  4  e1£))s(a1  4  a2  4  63) 

+  1/2  k1()  s\  4  1/2  kM  S2  4  1/2  kl(«3  -e,..;-  tj)2 
4  1/2  k2(63  4  e2  s  ax  -  62>2  4  1/2  kf  62 


+  2EX  f2  _  Vi  +  3  ^4 
e6  \  2  e6  e2 


(2.24) 


2ex 

where  g  is  the  gravitational  constant,  and  -  is  a  stiffness  factor  for 

e6 

the  beam.  The  strain  energy  of  the  beam  was  obtained  by  Integration  of 
the  square  of  second  derivative  of  the  deflection  curve,  from  Gq.  12.14). 
Similarly,  the  dissipation  energy  is. 


■  1/2  <^(63  -  6^  e^  c  -  ^)2  4  1/2  c2(63  4  6^  e2  c  -  $2) 
;2 


4  1/2  cr  63  . 


(2.25) 


2.3.3  Kinetic  Energy 

The  translational  kinetic  energy  is  obtained  from  velocity  of  the 
mass  center.  For  the  beam  element  an  Integration  of  the  square  of  velocity 
along  the  axis  of  beam  will  give  the  total  kinetic  energy.  As  outlined 
before,  this  energy  is  obtained  by  assuming  that  the  entire  mass  is  con¬ 
centrated  at  the  mass  center  of  the  beam.  The  angular  velocity  of  the 
mass  is  approximated  as 


V  s:  (&x  4  d2/2)K 


(2.26) 
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The  total  kinetic  energy  of  the  system  is  then  obtained  as 

T  -  1/2  m1  62  +  1/2  m2  fi2  +  1/2  MA  63  +  1/2  IA  a2 
+  1/2  ^[{63  s  ox  -  hx  «4  -  h2  d2  -  d1(eA  +  e,.)}2 

+  {63  c  a1  +  ^ (e3  -  hj^  6^  -  h2  c*2)}2J 

+  1/2  1^  +  d2/2)2 

+  1/2  Mct{53  a  -  6^  -  *1^4  +  eg)  +  +  ^2)  (e7  s(a2  +  ei^ 

-  eg  c(a2  + 

+  {«3  c  a1  +  “i<e3  -  $4)  “  (d^  +  d2)(e?  c^a2  +  91* 

+  eg  s(a2  +  ©x) ) >2 J 
+  1/2  +  d2)2 

+  1/2  Mp[{63  8  _  *4  "  ®l^e4  +  e6^  +  ^“l  +  ®2^e7  S^a2  +  91^ 

-  eg  c(a2  +  6^) 

+  «5  c(o2  +  6^  -  (d^  +  a2)(65  +  e15  -  eg  +  e10)s(o2  +  e^}2 

+  {«3  c  +  d^(e3  -  6^)  -  (dx  +  a2^e7  c^a2  +  91^ 

+  eg  s(a2  +  0^) 

+  d5  »(a2  +  0^)  +  (d^  +  d2)(«$5  +  e15  -  e^  +  e1Q) 

c(a2  +  e^}2] 

+  1/2  ID(d1  +  a2)2  . 


(2.27) 
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where 

1^  ■  moment  of  inertia  of  Body  A  along  K  direction  for  the  mass 
center  of  A,  etc. 

2.3.4  Generalized  Active  Force 

Now,  adopting  the  terminology  and  definition  of  Kane  [9],  one  needs 
to  find  generalized  active  forces  of  the  system.  Since  no  external  couple 
is  assumed,  forces  that  contribute  to  the  generalized  active  force  are 
the  breech  force  acting  on  the  gun  barrel  through  its  axis,  the  reaction 
force  from  the  recoil  mechanism  acting  on  the  housing  assembly,  and  the 
forces  f^  and  f2  acting  on  masses  and  n^,  respectively.  The  breech 
force  is  given  by 

F(t)«-F(t)  {c(a2  +  01)I1  +  s(a2  +  6^}  ,  (2.28) 

and  the  reaction  force  is  given  by 

R(t)  -  R(t)  {c(a2  +  9^  +  s(a2  +  01)31>  .  (2.29) 

Since  F(t)  can  be  considered  acting  at  the  mass  center  of  D,  one  has  the 
partial  rates  of  change  of  point  D*f  as 


R_D*  s 
v«  ■ 0 


W*.0 

®2 
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R-D*  7  ,  7 

V63  "  ®  «!  lx  +  c  ax  jx  , 


■  -  <«4  +  «65ll  +  (e3  -  54>5l 


-  eg  d-2  +  (6$  +  el5  -  eg  -  e?  +  e1Q)J2  , 


R-D*  7 

V*4"’il  ’ 


X*  ’  '  e8  *2  +  <S5  +  e15  -  e7  -  e9  +  e10>j2  • 


R-D*  7 

\mi2  * 


Hence,  contribution  of  the  breech  force  F(t)  to  the  generalized  active 
force  vector  is: 

[0,  0,  -  FCtJsCo^  +  cu  +  6^  , 

F(t){eg  +  (e^  +  ©g)c(c*2  +  6^  -  (eg  -  6^)a(a2  +  0^}  , 

F(t)  c(a2  +  ex)  ,  F(t)  eg  ,  -F(t)]T  .  (2.30) 


For  the  forces  f^  and  f2>  the  only  nonzero  components  are 


h 

F  A  -  f 

1  * 


F^  -  f 

\  *2 


(2.31) 
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respectively.  For  R(t) ,  the  partial  rates  of  change  of  positions  of 
and  Q2  are  the  same,  except  those  for  the  generalized  coordinate  6^,  and  a 
force  of  the  same  magnitude  acts  in  opposite  directions  at  Q^,  and  Q^, 
respectively.  Renee,  a  net  contribution  to  the  generalized  active  force  is 


F*(t)  -  R(t)  ,  (2.32) 

*5 

4  R-Q2  t  ,  R-Ql  t 

since  v^  •  i 2  and  v^  ■  0  . 


Therefore,  the  generalized  active  force  vector  of  the  system  in  R  is 
given  by 


-  -  F(t)  8(0^  +  <*2  +  0^  , 

Fa  -  F(t){eg  +  (e^  +  e6>c(a2  +  0^  -  (e3  -  64>s(o2  +  0^}  , 

F6  -  F(t)  c(ct 2  +  0X)  , 

4 

F  -  F(t)  eft  , 
a2  8 

¥t  -  R(t)  -  F(t)  .  (2.33) 


\ 

J  .  . .  .  -  . . . 
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2.3.5  The  Equations  of  Motion 

The  equations  of  motion  of  the  system  are  obtained  by  the  Lagrange's 
equation : 


d_  3T_  3(T  -  P)  3j£  p 

dt  3^  ‘  3q£  3qt  "  q± 


(2.34) 


where 


*1.  2, 


=  [6 


1»  o2*  °3*  “i*  °4*  a2 * 


V 


(2.35) 


As  seen  from  the  expressions  derived  before ,  the  kinetic  energy  of  the 
system  is  a  homogeneous  quadratic  function  of  the  variables  q^,  that  is. 


T 


1/2  H« 


(2.36) 


where  is  an  inertia  coefficient  of  the  system  in  R,  and  is  a  symmetric 
matrix  of  functions  of  the  generalized  coordinates  q^  only.  Therefore, 


(2.37) 


Defining 


Vv  V 


IF  -2? 
t  3^ 


3T 


(2.38) 


tA  repeated  subscript  indicates  simulation  having  a  range  of  1  to  7. 
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the  equations  of  motion  now  becooe 

H«  "j  +  ’j  ’  Gi  -  1'1-2 . 7  •  (2-39) 

Equation  (2.39)  can  be  put  in  a  first  order  form  as  follows: 

,  (2.40) 

Hij(qt)ai  *  °i(ui!'  V  ‘  Vv  •  (2-41) 

where  1,  J  -  1,  2 . 7. 

After  a  straightforward  differentiation,  the  elements  of  the  matrix 
H  and  vector  G  can  be  obtained.  Tables  1  and  2  give  expressions  for  the 
elements  of  G  and  H  in  the  original  notation  of  the  generalized  coordinates. 
For  abbreviation,  the  following  notation  is  introduced: 

h  -  e3  -  hx  -  h2  a2  , 

St  •  {5  +  *15  -  *9  +  *10  • 

e8  -  e7  s(a2  +  0^  -  eg  c(a2  +  0^  , 

ec  -  e7  c(a2  +  0^  +  eg  s(a2  +  0^  , 

et  -  e7  c^  +  o2  +  0j)  +  eg  s(a1  +  a2  +  0j)  , 
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-45  '  e4  +  e5  • 
e46  ‘  e4  +  e6  ’ 

V  V  ”l  •  hl  *4  -  h2  42  •  dl  e45  1 

v2  ”  ^3  c  “l  +  h  » 

v3  -  53  s  «i  -  54  -  e46  +  (a3  +  &2)es  , 

v4  ■  53  c  o3  +  -  5^)  ~  (a3  +  o2^ec  » 

v5  -  v3  +  55  c(a2  +  01)  -  (d3  +  a2)fit  8*°2  +  61*  * 

v6  "  v4  +  *5  S(a2  +  ®1)  +  (&1  +  d2)6t  C(“2  +  V  * 

To  complete  the  system  of  the  differential  equations,  initial  condi¬ 
tions  should  be  given.  As  noted  earlier,  the  weapon  is  assumed  to  be 
fired  automatically  at  the  equilibrium  configuration  under  gravity.  Hence 
the  initial  conditions  for  the  differential  equations  are  prescribed  such 
that  q1(0)  satisfy  the  equilibrium  conf iguratian  under  gravity,  and  the 
velocities  are  zero.  To  get  the  initial  configurations  one  has  to  solve 
the  static  problem 


3P 


0 


9 


i  m  1,  2,  . . . , 


7 


(2.42) 
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In  some  Instances,  one  may  Impose  certain  initial  values  to  specific 
generalized  coordinates;  for  example,  q^  *  6^  may  usually  be  specified 
initially  from  the  actual  measurement.  In  such  a  case,  the  number  of 
equations  in  Eq.  (2.42)  is  reduced,  accordingly.  The  equilibrium  configu¬ 
ration  will  be  denoted  by  a  superscript  e. 

3.  FORMULATION  AS  OPTIMAL  DESIGN  PROBLEM 


3.1  Introduction 

The  ultimate  goal  of  engineering  design  may  be  stated  as  that  of 
synthesizing  a  system  for  the  stated  needs  and  objectives.  A  conceptual 
model  for  the  system  is  usually  used  for  analysis  and  a  quantitative 
identification  of  the  objectives.  The  model  should  be  designed  in  such  a 
way  that  it  represents  the  system  as  closely  as  possible  and  is  still  with¬ 
in  the  analysis  capabilities.  In  the  previous  section,  such  a  model  was 
described  and  the  governing  equations  of  motion  were  derived.  In  this 
section  the  problem  is  posed  as  a  "parametric  optimal  design"  problem.  In 
this  problem,  the  maximum  of  cost  function  over  a  conceivable  range  of 
environmental  parameter  is  to  be  minimized  by  choosing  the  design  variables. 

3.2  Classification  of  Parameters 

The  first  step  in  an  optimal  design  algorithm  is  to  identify  various 
parameters  that  describe  the  system.  These  system  parameters  may  be 
divided  into  three  groups,  and  this  classification  is  described  for  the 
model  at  hand.  The  "state  parameters"  are  the  generalized  coordinates  of 
the  mechanical  system  and  these  are  6^,  f  ,,  6^,  a^,  6^ ,  a 2»  end  5^. 
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"System  design  parameters"  are  those  describing  the  system  geometries  and 
properties  of  the  specific  elements,  such  as,  k^,  k2Q,  k^,  k^,  c^,  C2> 

Cj. ,  El,  m^ ,  ^2 ,  M^,  ^ »  Ig >  J-g »  Ijji  *  ^2*  ®i*  e2*  on. 

Th.-  system  design  parameters,  to  be  determined  by  an  optimization 
technique,  are  chosen  by  the  designer.  The  last  group  of  parameters, 
called  the  "environmental  parameters,"  is  that  group  which  interacts  with 
the  environment  or  external  agencies.  Hence,  the  environmental  parameters 
can  be  characterized  as  having  some  degree  of  uncertainty  during  performance 
of  the  system.  For  the  present  model,  f^,  f^,  F(t)  and  R(t)  can  be 
considered  as  environmental  parameters.  It  may  be  noted  that  this  classi¬ 
fication  is  not  necessarily  unique,  but  depends  on  the  designer's  choice, 
especially  for  the  latter  two  classes  of  parameters.  For  example,  if  the 
recoil  mechanism  is  not  specified  explicitly  and  the  reaction  force  R(t) 
is  I.nown  to  vary  within  some  given  range,  one  could  choose  R(t)  as  an 
environmental  parameter  (in  this  case  one  may  call  it  an  environmental 
variable).  On  the  other  hand,  if  a  recoil  mechanism  is  to  be  designed  such 
that  the  system  performance  is  optimized  by  choosing  R(t) ,  then  R(t)  is 
classified  as  a  system  design  variable. 

In  the  following  formulation,  k^  (■  k2> ,  kf,  c^^  (-  Cj) ,  and  El  are 
chosen  as  design  variables,  while  the  other  system  design  parameters  are 
assumed  to  be  specified,  f^  and  f2>  related  to  the  road  conditions,  and 
related  to  operator's  choice,  are  chosen  as  the  environmental  para¬ 


meters 
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3.3  Environmental  Parameters  from  Road  Conditions 

It  is  assumed  that  the  road  shape  is  sinusoidal  with  amplitude  and 
wavelength  of  a  and  X,  respectively.  The  vertical  displacement  is  then 
given  by 


v  -  a  sin  “  t  (3.1) 

where  x  is  a  coordinate  in  the  direction  of  road.  Let  speed  s  of  the 
vehicle  be  uniform,  and  the  distance  between  front  and  rear  wheels  be 
b(“  e1  +  e2).  Then, 

x  -  st  ,  (3.2' 

vf  -  a  sin  yt  ,  (3.3) 

2tT8 

where  y  -  — —  ,  and  v  is  the  vertical  displ icement  of  the  rear  wheel 

iv  IT 

bottom.  Since  the  time  difference  between  the  rear  and  front  wheel  passing 
a  point  on  the  road  is  b/s,  the  vertical  displacement  of  front  wheel 
bottom  is 


v^  -  a  sin  y(t  -  b/s)  . 


(3.4) 


Denoting  the  time  difference  between  initial  firing  and  the  Instant  of 
zero  vertical  displacement  of  rear  wheel  bottom  by  f ,  one  obtains  expres¬ 
sions  for  the  vertical  displacements  (Fig.  2)  as 


1 
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a  sin  y(t  - 

f)  , 

(3.5) 

a  sin  y(t  - 

b/s  -  f)  . 

(3.6) 

It  say  be  noted  that  a,  y,  f,  and  s  are  the  environmental  parameters 
related  to  the  road  conditions.  Assuming  the  amplitude  is  also  fixed,  one 
has  y,  f,  and  s  (or  correspondingly  X,  f,  and  s)  as  the  environmental  para¬ 
meters. 

By  defining  and  such  that 

fl  *  k10vr  *  k10  a  8in  Y*fc  ~  f)  0*7) 

f2  *  k20Vf  “  k20  ®  8in  ~  b^8  "  *  (3.8) 

the  masses  «  and  t  .  of  the  model  on  an  even  surface  will  have  same 
response  as  that  of  the  vehicle  running  on  a  sinusoidal  road  [11]. 

3.4  Objective  Function 

As  described  earlier,  an  objective  for  the  system  must  be  formulated 
quantitatively,  possibly  as  a  function  of  the  parameters.  Sometimes,  the 
choice  of  a  function  form  for  the  objective  is  obvious,  but  in  general,  it 
depends  on  the  designer. 

For  the  weapon-vehicle  system,  it  is  assumed  that  the  motion  of  the 
system  is  induced  only  by  a  cyclic  breech  force  due  to  continuous  auto¬ 
matic  bursts,  starting  at  time,  t  ■  0.  Then  the  system  will  get  into  a 
steady  state  motion  after  a  transient  period.  Neglecting  the  factors 


such  as,  wind  velocity  and  gravitational  force,  trajectory  of  the  bullet 
can  be  assumed  to  be  a  straight  line.  A  plane  that  is  normal  to  the  axis 
of  gun  barrel  at  equilibrium  position  is  defined  as  a  target  plane.  A  set 
of  bullets  on  target  plane  at  a  distance  i  from  the  mass  center  of  gun 
barrel  is  shown  in  Fig.  3,  for  a  two-dimensional  case  and  for  a  one- 
dimensional  case.  For  general  environmental  conditions  with  Irregular 
road  shape,  the  shot  pattern  would  be  more  irregular.  In  any  case,  after 
firing  for  a  sufficiently  long  period  of  time,  the  shot  pattern  will  be 
bounded  by  a  finite  radius.  The  design  objective  is  to  have  an  optimal 
system  property  that  minimizes  the  disturbance  l.e. ,  the  radius  of  the  shot 
pattern  of  continuous  bursts.  The  idea  here  is  that  the  initial  few  shots 
are  more  Important  than  the  rest  of  the  shots.  Also,  the  deviation  of  the 
first  shot  from  the  aimed  point  is  not  important,  for  the  deviation  is  a 
constant  property  of  the  gun  and  can  be  assumed  fixed  for  a  given  gun. 
Hence,  without  loss  of  generality  one  can  assume  that  the  aimed  point  is 
the  same  as  the  point  where  first  bullet  lands. 

Now,  some  objective  functions  will  be  described.  The  choice  is  depen¬ 
dent  on  the  designer's  view.  The  simplest  one  may  be, 


J 

a 


N 


E-i*i 

N 


(3.9) 


where  r^  is  the  distance  from  the  aimed 
pattern,  and  is  a  weighting  factor, 
weighted  average  distance  of  the  first 
second  objective  function  may  be 


point  to  the  1th  point  of  shot 
This  objective  function  gives  the 
N  shots  from  the  aimed  point.  The 
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FIGURE  3.  SHOT  PATTERNS 
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FIGURE  4.  GEOMETRY  FOR  THE  TARGET  ANALYSIS 
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J 


d 


“i(ri+i 


» 


(3.10) 


denoting  the  weighted  sum  of  the  squares  of  the  consecutive  deviations  of 
the  first  N  shots.  This  form  of  the  objective  function  nay  be  good  If  a 
fast  steady  state  Is  needed.  A  third  objective  function  may  be 


i  - 


ra> 


(3.11) 


where 


r 

a 


r^N 


(3.12) 


Still  another  type  of  objective  nay  be  defined  as  a  functional. 


J 

c 


w(t)  r(t)dt 


(3.13) 


where  T  is  the  duration  of  operation,  u(t)  Is  a  weighting  function,  and 
r(t)  Is  the  distance  from  the  aimed  point  to  the  point  of  Intersection  of 
the  gun  barrel  axis  with  the  target  plane.  Other  objective  functions  are 
also  possible.  In  the  next  subsection  an  explicit  expression  for  r^  will 
be  derived  from  geometrical  considerations. 

3.5  Target  Analysis 

The  target  plane  was  defined  such  that  its  normal  was  along  the 
barrel  axis  at  equilibrium.  The  aimed  point  is  the  Intersection  of  the 
barrel  axis  with  the  target  plane.  In  Fig.  4,  the  following  vectors  are 
known,  after  the  equations  of  motion  have  been  solved: 
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fP/0  -  a*  I  -*•  a*  j  +  a*  K 

(3.14) 

rP'/0  -  ax  I  +  a2  3  +  a3  K 

(3.15) 

t6  e  -  e  -  e  - 

1  *  I  +  o2  J  +  otj  K 

(3.16) 

I  -  ax  I  +  a2  3  +  <*3  K 

(3.17) 

where  Ie  and  I  are  the  orientation  of  the  barrel  at  equilibrium  and  at 
arbitrary  time  t,  respectively,  and  (1,  J,  k)  are  same  as  defined  before. 
The  distance  l  of  the  target  plane  to  the  mass  center  of  the  gun  barrel  is 
known.  Then  from  the  geometry,  one  has 


_P'/P  _  jT/P  +  -T'/P  _  -T*/T  m  0  # 


(3.18) 


-T'/P'  T 
r  -  ai  , 


fT’/T  .  ie  .  0  . 


0.19) 


_T '  /T 

Let  the  unknown  vector  r  be  denoted  as 


fT  /T  ■  I  +  *2  J  +  X3  K  .  (3.20) 

By  substitution,  one  obtains  the  following  set  of  equations: 

"  ai  "  1  “i  +  '  "  0  > 


i  -  1,  2,  3 


(3.21) 
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Design  Variable,  b  =  (k^-  k2) ,  c^-  c2> ,  El]  ,  (3.27) 

Environmental  Parameter,  a  =  [8^,  a,  X,  f]  ,  (3.28) 

State  Variable,  z  £  [q^,  u^]  ,  m  -  1,  2,  3 . 7  (3.29) 

Then  the  optimal  design  problem  will  be  in  the  following  form: 

min  max  J  (3.30) 

oteA 

subject  to 

h(z,  b,  u)  -  0  ,  z(0)  given  ,  (3.31) 

and 

g(z,  b,  a)  <  for  all  ocA  ,  (3.32) 

where 

A  =  (a|q(a)  <  0}  .  (3.33) 

J  represents  one  of  the  objective  functions  from  Section  3. A,  and 
h(z,  b,  a)  "0  represents  the  equations  of  motion  from  Section  2.3.5. 
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This  formulation  is  more  general  than  the  conventional  nonlinear 
programming  problem  or  the  minimax  problem,  and  is  termed  as  "parametric 
optimal  design"  problem.  The  g-constraint  which  depends  on  the  environ¬ 
mental  parameter  a  is  called  "parametric  constraint."  The  a-constraints 
define  a  constrained  set  in  the  space  of  environmental  parameters.  Here 
one  may  choose. 


-  e  .  <  8,  <  0 

min  -  1  -  max 


s  < 


smax 


t 


X  .  <  X 

min  - 


> 


f2  <  (X/2s)2  , 


(3. 34) 

(3.35) 

(3.36) 

(3.37) 


where  the  subscripts  min  and  max  denote  the  given  lower  and  upper  bounds. 
The  last  inequality  indicates  that  the  time  difference  f  is  restricted  by 
the  period. 

A  transcription  of  the  problem  is  made  into  a  more  manageable  form. 

The  state  equations,  being  a  set  of  differential  equations,  are  discretized 
for  numerical  solution.  Also,  instead  of  the  min-max  type  of  objective, 
an  artificial  aasign  variable  bn+^  is  introduced.  A  transcribed  optimal 
design  problem  is 


”in  Vi 


(3.38) 
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subject  to 


max  JOE1,  b,  a)  <  b 
aeA 


n+1 


i+1  i  i  .  . 

z  -  z  ■  S(z  ,  b,  a) 


0 

z 


given  , 


(3.39) 


(3.40) 


g(b)  <  0 


(3.41) 


where 


z1  -  uo(ti)]T  »  °  "  1»  2 . 7  »  (3.42) 

and  t^  is  a  discretized  time  point  which  is  chosen  such  that  the  grid 
contains  the  firing  instant  where  the  values  of  z  are  denoted  by  an  upper 
bar.  S  is  symbolically  used  as  an  operator  giving  the  numerical  solu¬ 
tion  to  the  initial  problem.  In  Euler's  method,  this  would  be 


(3.43) 

Other  methods  of  solving  the  equations  of  motion  are  possible  and  can  be 
easily  implemented  in  the  solution  algorithm  for  "parametric  optimal 
design"  problem. 


S  =  At 


H"1  (G  -  Hu) 
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4.  SUMMARY 

In  this  report,  a  model  for  conceptual  design  of  an  automatic  weapon- 
vehicle  system  is  suggested  and  analyzed  such  that  a  study  of  the  inter¬ 
action  of  the  environmental  conditions  is  possible.  Although  the  motion 
is  constrained  in  a  plane,  the  model  has  several  important  features,  such 
as,  simulation  of  the  road  conditions  by  externally  varying  forces  f ^  and 
f2>  inclusion  of  the  control  of  the  azimuth  angle  which  may  be  randomly 
varied  by  the  operation,  simulation  of  the  recoil  mechanism  having  a  com¬ 
paratively  large  mass  moving  back  and  forth,  and  the  flexible  mount.  Due 
to  the  flexible  part,  the  model  has  infinite  degrees  of  freedom.  The  in¬ 
ternal  degrees  of  freedom  of  the  flexible  part  are  removed  by  the  assumption 
that  the  dynamic  configuration  is  the  same  as  the  static  configuration 
which  can  be  related  to  the  finite  number  of  generalized  coordinates  at 
the  interface.  Thus,  the  Infinite  degrees  of  freedom  system  reduces  to  a 
finite  degrees  of  freedom  system.  An  explicit  form  of  the  equations  of 
motion  are  then  derived  for  the  proposed  model.  The  system  has  seven 
generalized  coordinates  in  all,  and  the  governing  equations  of  motion  are 
highly  nonlinear  because  of  a  relatively  large  oscillating  mass.  The 
system  may  be  linearized  by  the  usual  method  of  linearization  in  the 
neighborhood  of  equilibrium  configuration  except  for  the  portion  involving 
large  oscillatory  mass. 

The  design  concept  adopted  is  to  find  the  design  variables  such  that 
the  system  has  a  desirable  property  of  least  possible  disturbance  from 
the  perturbation  due  to  a  continuous  burst  under  uncertainty  of  the 


* 
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environment  such  as,  road  conditions  and  the  azimuth  angle.  Some  expres¬ 
sions  for  this  objective  are  given  after  an  identification  of  the  para¬ 
meters.  In  the  formulation,  it  is  assumed  that  only  the  first  N  shots 
are  Important,  or  a  finite  duration  of  the  operation  from  the  start  of 
bursts  is  of  importance.  The  objective  functions  considered  are  functions 
of  only  the  position  of  the  mass  center  of  the  gun  barrel  and  its 
orientation. 

The  conceptual  optimal  design  problem  is,  then,  formulated  mathema¬ 
tically  as  a  "parametric  optimal  design"  problem.  A  numerically  implement- 
able  transcription  is  given  such  that  the  formulation  fits  into  the  finite 
dimensional  parametric  optimal  design  problem. 
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TABLE  1.  ELEMENTS  OF  R 


h12  -  ...  -  h1?  -  0  , 

H22  ‘  "2  ’ 

*  •••  *  ^7  ■  0  i 

H33-MA  +  MB  +  MC+MD  • 

H34  “  Vh  c  °1  "  *45  S  “l5 

+  (Mc  +  MD){(e3  -  6a)c  ax  -  c46  8  ^  -  *t> 

+  Mp  6fc  c(ax  +  a2  +  6^  , 

H35  '  -  «*1  -  Mc  '  V  “l  • 

tt36  '  '  h2  "b  0  »i  -  <»C  +  Vet  +  “d  {t  c(ol  +  *2  +  V  • 

“37  ’ "d  s(al  +  a2  +  V  • 

BW  ■  **  +  *B  +  h  +  h,  +  Vw  +  b2> 

+  (Mc  +  V< ««  +  <*3  *  V2  +  ‘7  +  4-  2\t  *.  -  2<*3  *  Vc1 

+  -  26t{e?  -  s(a2  +  0^  -  (e3  -  64)c(a2  +  0^}]  , 

H45  “  hl  ^  *45  +  vMC  +  **0^*46  *«^  +  5t  8  ^*2  +  61?  * 


38 


TABLE  2.  ELEMENTS  OF  G 


\ 

LJ 


\ 


G1  -  fx  +  cx(63  -  ej  c  b1  -  6^ 


‘  °1  8  "  k10  61  +  kl(63  ’  el  8  °1  "  61> 

G2  -  f2  +  c2(i3  +  ij,  o2  c  Oj  -  62) 

-  m2  g  -  k2Q  62  +  k2(63  +  e2  s  ax  -  62) 

G3  -  -  F(t)  6(0^  +  a2  +  0^  -  0^(63  -  c  -  i^) 

-  c 2(63  +  \  e2  c  ax  -  6J 

-  (MA  +  ”8  +  Mc  +  V8  -  kl(63  '  el  8  “l  -  V 

-  k2(63  +  e2  s  aj^  -  6^) 

G^  -  F(t){eg  +  e^fi  c(a2  -  0^  -  (e3  -  fiA>s(a2  +  0^} 

+  -  a1  e1  c  -  6j)c  0^ 

-  c2  e2^3  +  e2  c  al  “  °1 

+  Mg  ^(Vj^  c  <Xj  ~  v2  8  ax) 

t 

+  Mc  53(v3  c  8  a j) 

+  Mg  «3(v5  c  ox  -  vg  s  ax) 

-  Mg  g(h  c  al  -  eA5  s  c^) 

-  (Mc  +  MD)g{(e3  -  64)c  ax  -  eA6  8  ^  -  et> 

"  *D  8  6t  c(al  +  a2  +  91) 

+  kj^  -  e1  8  ^  -  «1)e1  c  0^ 

-  k2(63  +  e2  c  a1  -  «2)«2  c  0^ 


»  / 

Precedlni  page  hlank 
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TABLE  2  (cont'd) 


G5  -  F(t)c(a2  +  0^-013  hx  v2  +  Mc  v4  4  Mp 

+  hl  ^  §  8  al  +  (\  +  Mpjg  8  ai 


.  6EI  12EI  - 

+  —  °2  — r  6  a 
e6  e6 


G6  "  F(t)e8  "  “b  V2  h2  “l  +  MC(“l  +  d2)<V3  *c  +  VA 

+  +  a2){v5  ec  ~  6t  V5  c(a2  +  6l>  +  v6  es 

-  Mp  «5{v5  s(a2  4  0^  -  vfi  c(a2  4  0^} 

4  M„  g  h„  s  a,  4  (M„  4  M_)e 

D  A  X  U  U  t 

“  «jj  g  0(0^  4  a2  4  01)5t 


AEI  .  6EI  . 
~t  “2  +  5« 
6  *6 


?7  »  R(t)  -  F(t)  -  cr  «5 

4  Mp  4  d2){c(a2  4  0^)  -  s(a2  4  0^)) 

-  Mp  g  s(a1  4  a2  4  0^  -  kf  65  . 


-  5t  v6  s(“2  +  V1 


